Stationary point approach to the phase transition 
of the classical XY chain with power-law interactions 
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The stationary points of the potential energy function V of the classical XY chain with power law 
pair interactions (i.e., interactions decaying like r~ a with the distance) are analyzed. For a class 
of "spinwave-type" stationary points, the asymptotic behavior of the Hessian determinant of V is 
computed analytically in the limit of large system size. The computation is based on the Toeplitz 
property of the Hessian and makes use of a Szego-type theorem. The results serve to illustrate 
a recently discovered relation between phase transitions and the properties of stationary points of 
classical many-body potentials. In agreement with this relation, the exact phase transition potential 
energy of the model can be read off from the behavior of the Hessian determinant for exponents 
q between zero and one. For a between one and two, the phase transition is not manifest in the 
behavior of the determinant, and it might be necessary to consider larger classes of stationary points. 
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I. INTRODUCTION 

As has long been known, the stationary points of a clas- 
sical Hamiltonian function or potential energy function 
can be employed to calculate or estimate certain physi- 
cal quantities of interest. Well-known examples include 
transition state theory [l[ or Kramers's reaction rate the- 
ory for the thermally activated escape from metastable 
states [2| , where the barrier height (corresponding to the 
difference between potential energies at certain station- 
ary points of the potential energy function) plays an es- 
sential role. More recently, the noise-free escape from 
quasi-stationary states (i. e., metastable states whose life- 
times diverge with the system size) has been related to 
the presence of stationary points of marginal stability 
0. Apart from studies of dynamical properties, station- 
ary points have also been extensively used for estimating 
thermodynamical properties by means of the superposi- 
tion approach 

Dynamical properties like the aforementioned ones are, 
as one might expect, not unrelated to the statistical phys- 
ical behavior of a system. Accordingly, as worked out 
beautifully in Q, properties of stationary points of the 
potential energy function Q reflect in dynamical and sta- 
tistical physical quantities simultaneously. This observa- 
tion sparked quite some research activity, reviewed in Q , 
with the aim of relating equilibrium phase transitions to 
stationary points and their indices. Most importantly, it 
was shown in Q that, under a number of technical con- 
ditions, the presence of stationary points of the potential 
energy function is necessary for a phase transition to take 
place. Subsequently, it was noticed in |9( that the Hessian 
determinant of the potential energy function, evaluated 
at the stationary points, adds a crucial piece of informa- 
tion for discriminating whether or not a phase transition 
occurs. Omitting some of the technicalities, the essence 



of the criterion on the Hessian determinant can be cap- 
tured as follows [Toj | . 



Criterion. Let 



C 



N 



H (P,V) = jJ2pl + V{q 1: ...,q N ), 



(1) 



k=l 



with some constant C 0, be the total energy func- 
tion of a system with N degrees of freedom, where p — 
(pi, . . . ,Pn) and q — (qi, . . . , <7jv) denote the vectors of 
momenta and positions. The potential energy V will, in 
general, have stationary points q^ defined as solutions of 
the set of equations 



dV(q) 



dqk 



k = 1, 



,N. 



(2) 



The stationary points are assumed to be isolated, and 
their number is assumed to grow at most exponentially 
with N. In the thermodynamic limit N — > oo, the sta- 
tionary points can induce a phase transition at some crit- 
ical energy per degree of freedom e c only if the following 
two conditions are met: 



1. There exists a sequence {<7;f }^ =JV 
points of V such that 



lim 

AT->oc 



N 



of stationary 



(3) 



converges and v c — (v)(e c ) is the ensemble expec- 
tation value of v = V/N at the energy e c . 

2. The asymptotic behavior of the Hessian matrix ri 
of V , evaluated at the critical points q^ contained 
in that sequence, is such that 
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lim \det rl(q 



= 0. 



(4) 
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In short, the criterion requires the existence of a se- 
quence of stationary points whose potential energy con- 
verges to v c and whose Hessian determinant vanishes in 
the sense of (UJ) in the thermodynamic limit. Setting the 
constant C in (JlJ to zero, the kinetic energy term is ab- 
sent, as is the case for many classical spin models. The 
criterion therefore remains valid in this case, with the 
only differences that the total energy H equals the po- 
tential energy V, and the critical energy e c is identical to 
the critical potential energy v c [TTI | . This is the case we 
will be concerned with in the present article. 

Note that the above criterion is not sufficient for a 
phase transition to occur: finding a sequence of station- 
ary points with the behavior specified above does not 
guarantee a transition to take place at the corresponding 
critical energy. However, as model calculations suggest, 
the criterion usually appears to single out the correct 
transition energies [9j. Importantly for the application 
of the criterion, it is not necessary to know all stationary 
points of V, but a suitably chosen subset may be suffi- 
cient. This matter of fact was pointed out by Nardini and 
Casetti, and suitably constructed sequences of stationary 
points were used in [To| to single out the phase transi- 
tion of a model of gravitating masses and analytically 
determine its critical energy. 

Comparing the Hessian determinant criterion to other 
analytic tools in the statistical physics of phase transi- 
tions, its remarkable property is that it is local in con- 
figuration space. In contrast to, say, the calculation of 
a partition function, no averaging over a large, high- 
dimensional manifold is necessary. Instead, only the lo- 
cal properties of a sequence of stationary points need to 
be analyzed. Of course, finding an appropriate sequence 
of stationary points can be equally hard or impossible, 
but in certain instances such a local approach may prove 
beneficial. 

In this article, we study the stationary points of the po- 
tential energy function of a chain of classical XY spins 
(or rotators), coupled by a pair interaction which decays 
like r~ a with the distance r of the spins on the lattice. 
The model is introduced in detail in Sec. HI1 Although 
one-dimensional, it shows a phase transition from a fer- 
romagnetic to a paramagnetic phase for exponents a be- 
tween zero and two, and the aim of the present work is to 
explore the relation between stationary points and phase 
transitions for these values of a. 

There are a number of interesting aspects of this study 
that deserve mention: First, inspired by Ref. (io| . in Sec. 
Illl a method is devised of how to construct special classes 
of stationary points for lattice spin systems. The poten- 
tial energy at such stationary points is evaluated in Sec. 
IVI The Hessian at such a stationary point, as required 
by is found to be a Toeplitz matrix. As carried out 
in Sec. El this property allows us to employ a Szego- 
type theorem for the calculation of the asymptotic be- 
havior of the Hessian determinant in the limit of large 
system sizes N. The results of Sec. fVl depend on the ex- 
ponent a not only quantitatively, but also qualitatively: 



For ^ a ^ 1, the asymptotic behavior of the Hessian 
determinant indeed signals the phase transition at the 
exact value of the transition energy, as purported by the 
criterion of Sec. HI For I < a ^ 2, no signature of the 
phase transition is detected from the Hessian determi- 
nant of the special class of stationary points considered, 
and one might conclude that other (or larger) classes of 
stationary points have to be taken into account. The 
findings are summarized and discussed in Sec. IVI1 



II. CLASSICAL XY CHAIN WITH 
POWER-LAW INTERACTIONS 

Consider a set of N lattice sites labeled by an inte- 
ger number j 6 {I, . . . , N} where, to ease the notation, 
we assume N to be odd. To each site, a planar vector 
of unit length is assigned, parametrized by the angular 
variable 6>j 6 (— ir, tt]. The classical XY" chain with power 
law interactions is characterized by the potential energy 
function 



jV (N-l)/2 

i=i j=i 



1 - cos(6»i - 9 i+ j) 



(5) 



with 9 = (#i, . . . , 9pf) and some nonnegative exponent a. 
Although suppressed in the notation, indices i of the 6i 
variables are always to be considered modulo N, such as 
to account for periodic boundary conditions and to guar- 
antee indices in the range from I to N. The potential 
energy function ([5]) describes N classical spin variables 
on a ring (chain with periodic boundary conditions), 
where each spin interacts with every other. The inter- 
action strength between two spins decays proportionally 
to l/j a , where j is the minimal distance of the two spins 
on the ring. The potential energy (JSJ) is endowed with a 
normalization factor defined as 



(6) 



The asymptotic behavior of Af in the limit N — > oo can 
be computed, yielding 




(I - a)2 1 - Q A^"- 1 
2.V - I 1/lniV 
I/C(a) 



for < a < 1 , 
for a = 1 , 
for a > I, 



(7) 



where C denotes the Riemann zeta function. This nor- 
malization factor, introduced in Ref. [12], is chosen such 
as to guarantee extensivity of the potential energy, i.e., 
a finite limit of the potential energy per particle in the 
limit N — > oo. 

The thermodynamic behavior of this model depends 
on the exponent a in the following way: For ^ a ^ 
1, the thermodynamic behavior is identical to that of 
the mean-field (or Curie- Weiss) case a — 0, showing a 



ferromagnetic continuous phase transition characterized 
by mean-field critical exponents 13]. For 1 < a ^ 2, the 
model also shows a phase transition, but thermodynamic 
functions differ from the mean-field case (see [1J] and 
the comment in Sec. 5 of Ref. [HI). For a > 2, no phase 
transition occurs. The three regimes for the exponent a, 
the corresponding thermodynamic behavior, and also the 
methods of proof are analogous to Dyson's analysis of the 
Ising chain with spin-spin interaction strengths decaying 
as l/j a [3. 

III. STATIONARY POINTS 
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FIG. 1. Sketch of stationary points of the potential V for 
N = 8, where is the angle between the arrow and the dashed 
axis. Top: Stationary points where all Oi £ {0,7r}. Bottom: 
Spinwave stationary point @, where all differences Oi — 6i-i 
between neighboring angles are equal, with differences chosen 
such that Go = On, in compliance with the periodic boundary 
conditions. 



Stationary points of the potential energy function ([5]) 
are defined as the real solutions of the set of equations 







dV(9) 



E 

3=1 



r 



for k — l,...,N. Since the potential energy function 
([5]) is invariant under a global rotation Oi — > 9{ + <f> with 
4> £ R, the solutions of ([H} come in one-parameter fam- 
ilies: Given a stationary point (#1, . . . , On), every point 
(61 +</>,..., On + <p) is also a solution of ©. The cri- 
terion of Sec. HI however, requires all stationary points 
to be isolated, and we therefore have to get rid of the 
trivial rotational degeneracy. We explicitely destroy the 
rotational symmetry by fixing On — and eliminating 
the equation with fc = N in (|8[). The thermodynam- 
ics of this reduced model is identical to that of the full 
one, since the contribution of one degree of freedom to 
the partition function is negligible in the thermodynamic 
limit. 

Determining all solutions of the remaining set of non- 
linear equations is presumably a hard task (too hard for 
the author at least). There are, however, two particu- 
larly simple classes of solutions, similar in spirit to those 
constructed in [lol ] for a one-dimensional model of grav- 
itating masses: First, any combination of Oi £ {0, n} for 
i = 1, . . . , N — 1 will make the sine functions in ([8]) van- 
ish and therefore satisfies the set of equations. A second 
class of solutions is given by 



Q$ = mx with x — 2im/N, 



(9) 



where m, n £ {1, . . . , N} and hence < x < 2tt. These 
solutions have constant radian x between neighboring 
spins, implying sin(#fc — Qk+j) — sin(#k_j — Ok) and there- 
fore each of the summands in ([8]) vanishes separately. 

These two classes of solutions are illustrated in Fig. Q] 
and, as is easily checked numerically, are not exhaustive. 
This is probably expected, in particular when comparing 
with the results for the nearest-neighbor XY chain for 
which all stationary points can be computed analytically 
1171. The two classes of solutions introduced above are 



also solutions in the case of nearest-neighbor interactions, 
but many more solutions exist. Since nearest- neighbor 
interactions can be considered as the limit a —¥ 00 of 
the power law decay discussed in the present article, it is 
maybe not too surprising to find that (at least many of) 
these solutions persist to finite a. 

We will in the following restrict the analysis to "spin- 
wave" stationary points ©, mainly for the reason that 
the Hessian matrix at these points, as discussed in de- 
tail in Sec. |V) is a Toeplitz matrix. This structure is 
particularly helpful when calculating the large- N asymp- 
totics of the Hessian determinant. Moreover, from the 
results on the nearest-neighbor XY chain in pj}, one 
may be led to conjecture that the spinwave stationary 
points arc of particular importance for our purposes: At 
least for nearest-neighbor interactions, the spinwave sta- 
tionary points have the smallest absolute value of the 
Hessian determinant amongst the stationary points of a 
given potential energy, and therefore determine whether 
the criterion of Sec. U is satisfied or not. 



IV. POTENTIAL ENERGY AT STATIONARY 
POINTS 

The criterion of Sec. Q] involves the potential energy 
evaluated at the stationary points. Inserting the spin- 
wave stationary points (0) into the potential energy func- 
tion ([5]), we obtain 



v(x) 



V(0 (x) ) N 



N 



N 



N (N-l)/2 

E E 

i=l j=\ 



1 — cos(jx) 



(N-l)/2 



cos(jx) 



3=1 



Since I cosed € 1 for all x, we have 



(N-l)/2 

« E 

3=1 



cos(jx) 



r 



(JV"l)/2 , 1 

^ E i4> 



3 = 1 



3 



(10) 



(11) 



confirming that the normalization factor j\f in © had 
been chosen appropriately in order to render the poten- 
tial energy per spin finite in the thermodynamic limit. 
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In the limit N — > oo and for certain values of the expo- 
nent a, the summation in the second line of (|10l) can be 
performed explicitely: For a = 1 we use formula 1.441.2 
of OH, and for a e 2W formula 1.443.1 of the same 
reference, to obtain 



1 JV Jln[2(l-cosx)] 
V{X) - 2 + T i (-47r 2 )«/ 2 B Q [x/(27r)]/a! 



for a = 1, 
for a e 2JN , 

for the potential energy of a spinwave solution in the 
thermodynamic limit. B a denotes the Bernoulli polyno- 
mial of order a as defined, for example, in Sec. 9.62 of 
[Hjj. The graph of v(x) is shown in Fig. [2] (upper plot) for 
exponents a — 0,1,2,4, and 6. For noninteger values of 
a, the infinite sum in (1101) cannot be performed, but v(x) 
can be evaluated numerically for reasonably large system 
sizes N. The resulting curves (not shown in Fig. [2]) are 
found to interpolate smoothly between the curves for in- 
teger a. For a given positive, even a = 2, 4, 6, ... , the 
potential energy values cover densely the entire range of 
potential energies per spin accessible to the system in the 
thermodynamic limit. This is a desirable property when 
applying the criterion of Sec.Hl as it allows us to use spin- 
wave stationary points for the construction of sequences 
of stationary points whose potential energies converge to 
any desired value accessible to the system. 

For exponents in the range ^ a ^ 1, the situation is 
more intricate. For large, but finite, system sizes N, the 
potential energies corresponding to the spinwave station- 
ary points Q^ x > become denser and denser on the interval 
(0, 1/2) with increasing N. In the thermodynamic limit, 
however, the potential energy converges to 1/2 for any 
given value of x G (0,27r), and to zero for x — 0, re- 
sulting in the straight line plotted in Fig. [2] (upper plot). 
The approach to this behavior with increasing system 
size N is illustrated in Fig. [5] (lower plot) for the expo- 
nent a = 1/2. 

Equation ([12"]) and the numerical results in Fig. [2] 
demonstrate that the potential energy per spin of a se- 
quence {8^} of spinwave stationary points for increas- 
ing N, but with a fixed value of x, indeed converges to 
a limiting value, as required in equation (J3j) of the cri- 
terion in Sec. HI To construct a sequence of stationary 
points with a given value of x, it will in general be nec- 
essary to restrict the sequence to some infinite subset 
{Ni,N2, ■ ■ ■} of system sizes such that 2im/Ni — x for 
some n G {1, ... , Ni}. 



at a stationary point. Again, in order to destroy the 
trivial global rotational symmetry of the potential energy 
function V, one of the spin variables, say, 9n, is fixed 
at zero. The resulting potential energy is a function of 
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FIG. 2. (Color online) Upper plot: The graph of the poten- 
tial energy per spin e(x) of a spinwave solution 9^ x > in the 
thermodynamic limit as given in (|12[) . For a given positive, 
even a = 2, 4, 6, ... , the energy values cover densely the en- 
tire range of energies per spin accessible to the system in the 
thermodynamic limit. Lower plot: Potential energy per spin 
e(x) for a — 1/2 and spinwave stationary points |9|, plot- 
ted for various system sizes N. With increasing N, the curve 
approaches a horizontal line of energy 1/2. 



V. HESSIAN DETERMINANT AT 
STATIONARY POINTS 



In order to apply the criterion stated in Sec. U we also 
need to evaluate the determinant of the Hessian matrix 



N — 1 variables 6±, . . . ,9n-i, and its Hessian Hn is an 
(N — 1) x (N — 1) symmetric matrix with entries 



5 



[Hn]h(0) 



d 2 V{6) 
d6 k d6i 



(jV 12 CQS (-g fc _ + cos (g fc _ ^ ^ 



4=1 

COs(flfc - flj) 

A(Z- fc) Q 



(13) 



for k I, 



for k, I = 1, . . . , N — 1, where 



A(Z-fc) = 



for |i-jfe| (AT-l)/2, 



2V - |Z - ftl else, 



(14) 



is the minimal distance between k and I on the ring. 
Evaluating the Hessian at a spinwave stationary point 
0<» as defined in ©, one obtains 



1 - 2t>(z) 



for k — I, 
for k ^ I, 



with 



,W 



-A/" 



cos(jx) 
A(j) Q 



(15) 



(16) 



and with the potential energy per spin v(x) as given in 
(ITU)) . Without fixing 6^ to zero, this matrix would be 
circulant and the eigenvalues were readily obtained by 
Fourier transforming a row vector of the matrix. Fixing 
On corresponds to eliminating the TVth row and column 
of the matrix, and although the resulting matrix is not 
circulant anymore, it retains the Toeplitz property: As 
is evident from (|T51) . the elements [H.N]ki depend only on 
the difference I — k of the indices. 



A. Szego's theorem 

For our purposes, the Toeplitz property comes in 
handy, as a number of theorems on the large- N asymp- 
totics of determinants are known for sequences of N x N 
Toeplitz matrices [lj|. The kind of sequence {Ijv}jv =1 
of matrices Tjv that is typically considered in the math- 
ematics literature is where the matrix elements 



[T N (f)} kl =U_ k (f) 



(17) 



are given as Fourier coefficients of a complex- valued func- 
tion / defined on the circle, 



1 

27 



For particularly well-behaved /, Szego's theorem states 
that the large- N asymptotic behavior of the determinants 
of such a sequence is given by 

Jim |det(2V(/))| 1/JV - exp (A. / hxf^dA , 

(19) 



and many generalizations of this result to larger classes 
of symbols / can be found in the literature [19( . Inverting 
the Fourier transformation (118[) . we can write 



(20) 



j = -oo 



For the XY chain, using equations (|T5j) and (jlpj) and a 
standard trigonometric identity, we obtain 



fW(4>) =l-2v(x) 



oo 

4=1 



cos[j(:r + cj))} + cos\j(x - 0)] (21) 



for the symbol of the Hessian, evaluated at a spinwave 
stationary point 0< x ) as defined in ©. Then, as in the 
calculation of the potential energy in Sec. lIVl the formulae 
1.441.2 and 1.443.1 of [18j can be used to perform the 
summation in (|2"Tj) for the values a = 1 or a G 2 IN of the 
exponent. 



1. Exponent a — 1 



In the case of a = 1 we can use the identity 



E"° cos(jx) 1 , .. 
—f-L = --h 1 [2(l-COBx)] 



4=1 



(formula 1.441.2 of [18]) to write jUJ in the form 



/<*>(«£) = y In 



cos x — COS < 
COS x — 1 



(22) 



(23) 



To compute the Hessian determinant as a function of 
the potential energy per spin, we invert the first case in 
equation (fT2|) . yielding 



r\ 1 1 ( 2v ~ 1 

cos x(v) = 1 exp — —. — 



(24) 



(lg) Inserting this expression into (|23p gives 



/(*(«))(<£) = l-2w+AAln 



2(1 — cos (j>) — exp 



2v-l 



(25) 

As a consequence of the asymptotic behavior ([7]) of TV, 
the logarithmic term on the right hand side vanishes and 
we obtain 



/(*(«)) (0) =l-2u. 



(26) 
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\detH^ 1/N 




plotted in Fig. [3] together with numerical results for the 
Hessian determinants for several finite system sizes. The 
convergence of the finite-system data to their infinite- 
system limit T>i is slow, but this is no surprise as the 
finite- N corrections in (|25|) are logarithmic. 

Interpreting (|2"T1) in terms of the criterion of Sec. Q] 
[and equation ((¥]) in particular], we observe a vanish- 
ing Hessian determinant at the potential energy per spin 
v = 1/2. This value coincides precisely with the known 
phase transition potential energy of the model and nicely 
illustrates the criterion of Sec. fl] 



FIG. 3. (Color online) The iVth root of the Hessian deter- 
minant det"Hiv for a = 1, evaluated at spin wave station- 
ary points 6^ x \ plotted versus the corresponding potential 
energy per spin v(x). With increasing system size N, the 
numerically computed determinant slowly (logarithmically) 
approaches the analytic large- N asymptotic result T>\ (black 
line). 



From equation (TT9)) , the large- N asymptotic behavior of 
the Hessian determinant (TT51) is found to be 



X>i(u) := lim I det (6 ,(a:) ) 1 1/A " = l-2v (27) 



in the case of a = 1, valid for accessible potential energy 
values from the interval [0, 1/2]. The straight line (f2~Tf is 



2. Exponents a G 2JN 
In the case a £ 2Mo we can use the identity 

£=^~c-i>-i^.(£) w 

3 = 1 J 

(formula 1.443.1 of pj). valid for s$ x ^ 2tt. From this 
formula, and considering that x, 4> € [0, 27r], we can write 



E 



cos\j(x±4>)] 



_ ( _ ir / 2 l(2 7 r) a 



B n 



x±4 

2tt 



2 al 
for ^ x ± d> < 2tt, 



< > ' \ (29) 
^y-^rTl) else. 

Inserting these identities into (|2ip. we obtain 



,7V (27T)" 



B a 




B a \ 


x+4> 

2ir 


B a \ 


x+tp 
2tt 


B a \ 


x+4> 
2tt 



-2B a {^ ) 



2tt 



lj -2B a (£)+B a 
— 2B a (^) + B a 



x — <p 
2tt 



X — q> 

2tt 



for (p ^ x and (j) ^ 2ir — x, 
for <fi ^ x and <j) > 2tt — x, 
for (j) > x and <j) > 2~k — x, 
for d> > x and 6 ^ 2ir — x. 



(30) 



For our purposes, the case a = 2 is particularly interesting, as this is the only positive even exponent for which the 
XY chain with power law interactions exhibits a phase transition. In this case, making use of the Bernoulli polynomial 



B 2 (x) 



1/6, the symbol simplifies to 



-2n) 2 
+ 2n(x - 



for d> < x and ch < 2ir 



x, 



2ir) 2 — 2ir(x — (f)) for <p ^ % and (j) > 2ix — x, 



for (f> > x and <j) > 2~k — x, 
for d> > x and ci ^ 2tt — x. 



(31) 



To compute the Hessian determinant as a function of the 
potential energy per spin, we invert the function v(x) in 
([i~2|l . yielding 

x(v) = 7T (l ± y/l-4v/3) (32) 



in the case of a — 2, and insert this expression into (|31l) . 

The symbol (</>), plotted in Fig. HI is easily seen 
to be positive for some values of x and </>, and negative 
for others. In principle this causes a problem when com- 
puting the large- N asymptotic behavior of the Hessian 
determinant (1151) from Szego's theorem (|19|) where we 
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FIG. 4. (Color online) The graph of the symbol f (x) (<i>) for 
a — 2 as given in (|31[) . plotted as a function of x and <f>. 



have to integrate the logarithm of f^. Instead, to cir 
cumvent this problem, we chose to replace in (pjj 
by its absolute value and compute 



V 2 := exp 



(^J^H/ (x( "%)rV). (33) 



Inserting (f3"Tj) and performing the integration, we obtain 



^2 = -(3-V9-12w)' 



exp 
exp 



-2 + 2J2y/l -4v/3- larccoth 



3-4-u 



2 + 2J2-y/l - Av/3 - 1 arctanh ( 



2V9-12U-3 



3-4-u 
2V9-12t>-3 



for v < 9/16, 
for 9/16 < t> < 3/4. 



(34) 



IdetWjvl 1 ^ 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 



FIG. 5. (Color online) The Afth root of the Hessian deter- 
minant detHN for a — 2, evaluated at spin wave stationary 
points 6^ x \ plotted versus the corresponding potential en- 
ergy per spin v(x). Already for moderate system sizes N, the 
finite-system data are in excellent agreement with the analytic 
large- A^ asymptotic result T>2 (black line). 



The graph of T>2 is shown in Fig. [5] together with nu- 
merical results for the Hessian determinants for several 
finite system sizes. The numerical results are in such ex- 
cellent agreement with (|3~4")l that it is tempting to believe 
that taking the absolute value of f( x ( v » in (|3"3"1) is not 
merely an approximation, but gives an exact asymptotic 
expression for the Hessian determinant. Unfortunately, 
the author was unable to proof this conjecture (2pj . 

Interpreting (|3~4"|) in terms of the criterion (Q| , we ob- 
serve a strictly positive Hessian determinant on the entire 
range of accessible potential energies per spin v € [0, 3/4]. 
For the class of spinwave stationary points considered, 
the analysis of the Hessian determinant therefore fails to 
give an indication of the phase transition known to exist 



for a — 2. 

For even exponents a = 4, 6, 8, ... , analogous calcu- 
lations can be performed. The corresponding asymptotic 
results for the determinant, computed according to (|33p . 
share the most important features of the case a = 2: The 
analytic large-A^ asymptotic result T> a is bounded away 
from zero (see Fig. [5]) and in excellent agreement with 
numerical data (not shown). For a > 2, however, the 
fact that T> a is bounded away from zero was to be ex- 
pected, as no phase transition occurs in this case. The 
graphs in Fig. [5] also suggest that, in the limit a — > oo, 
T> a approaches the function X>oo(^) = \2v — 1|. This limit 
corresponds to nearest-neighbor interactions on the lat- 
tice, and indeed coincides with the behavior of the 
determinant of the nearest-neighbor XY chain reported 
in Ref. [17] . Despite the absence of a phase transition 
in the XY chain with nearest-neighbor interactions, Poo 
vanishes at v = 1/2. Note that this finding is not in con- 
flict with the criterion of Sec. HI since a vanishing D is 
not claimed to be sufficient for a phase transition. More- 
over, as explained in more detail in Ref. [l7T ]. the poten- 
tial energy v = 1/2 at which vanishes corresponds 
to infinite temperatures in the nearest-neighbor model, 
and at least within the standard canonical setting where 
temperature is the control parameter, such a transition 
would be elusive in any case. 



3. Other values of a 

For other values of the exponent a, Szego's theorem 
can still be used to obtain the asymptotic behavior T> a of 
the Hessian determinant but, to the best of the author's 
knowledge, the infinite sums in (fTU|) and (|2ip cannot be 
calculated anymore. Numerical results for the cases a = 
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FIG. 6. Asymptotic large- JV behavior T> a of the Hessian de- 
terminant, plotted as a function of the potential energy per 
spin v for various exponents a. Dark lines correspond to val- 
ues of a for which exact analytic results are available, light 
lines (for a = 6/5 and a = 3) were obtained numerically. 
The dark solid straight line is obtained analytically for a = 1, 
but coincides with the numerical results [and also with the 
Hadamard bound (|39|) ] for a £ [0, 1]. T> a vanishes at v = 1/2 
for a £ [0, 1] (long-range interactions) as well as for a = oo 
(nearest-neighbor interactions). For all finite a > 1, T> a is 
bounded away from zero. 



5/4 and a = 3 are shown in Fig.[6l Again, T> a is bounded 
away from zero for all a > 1. 



Hadamard bounds 



The Szego-type theorem we used in Sec. IV Al allowed us 
to obtain exact asymptotic large- TV results of the Hessian 
determinant of spinwave stationary points. The draw- 
back, however, is that an evaluation of the resulting infi- 
nite Fourier sums is possible only for the exponents a = 1 
and a £ 21No. In the present section we will supplement 
these results by an upper bound on T> a , valid for any 
a ^ 0. Comparing to numerical data, we will observe 
that the bound is sharp for ^ a ^ 1. 

A bound on the determinant of a (real or complex) 
N x Af-matrix M can be obtained by the celebrated 
Hadamard inequality 



N 



| det Af | < J] 



(35) 



where Cj denotes the jth column (or row) vector of M , 
and ||cj|| its Euclidean norm. In contrast to the methods 
of Sec. IV Al that are based on the Toeplitz property of the 
Hessian, the Hadamard inequality can be used to bound 
the Hessian determinant not only of spinwave stationary 
points ([9]), but also of any kind of stationary point. In 
the context of phase transitions and their relation to sta- 
tionary points and their determinants, Hadamard bounds 
first have been used in [lOj . 

The Hadamard bound (|35|) becomes particularly sim- 
ple for a circulant matrix. In this case, ||c,-|| = ||cfc|| for 



all j, k = 1, . . . , JV, and hence 

\detM\^ N s: ||cj|| 



(36) 



for any j. The Hessian (II 5[) of a spinwave stationary point 
we want to study is not quite a circulant matrix, but it is 
closely related: Hn(0^ is an (JV — 1) x (JV — l)-matrix, 
obtained from an JV x JV circulant matrix by deleting 
one row and one column. The norm of every column of 
Hn(9^ is therefore bounded above by the norm 



JV 



INI 



J=2 



os 2 [x{j - 1)] 

A (j - 1)2« 



1/2 



(37) 



s$ y/[l-2v(x)] 2 +M 



of any column of the original JV x JV circulant matrix. A 
bound on the Hessian determinant is then given by 

|detftiv(0 (a:W) )| 1/JV < ||c|| < v/(l-2v) 2 +AA. (38) 

For a £ [0, 1], we have observed in ([7]) that the normal- 
ization constant Af goes to zero in the limit JV — > oo, 
yielding 



lim \detH N (0 ix(v)) )\ 1/N < \l-2v\. 



(39) 



For the accessible values of the potential energy per spin 
v £ [0,1/2], this bound coincides with the exact a = 1 
asymptotic result for T>i obtained in Sec. IV Al Compar- 
ing (|39[) to numerical data, it is tempting to conjecture 
that the bound is tight, i.e., coinciding with the exact 
asymptotics, for all a £ [0,1], but the author was not 
able to prove this. Such a result for T> a (i. e., one that is 
independent of the precise value of a) would also agree 
well with the known fact that the thermodynamics of the 
XY chain has no a-dependence as long as a is between 
zero and one [l3j]. In particular, the bound (|3"9")l vanishes 
at v = 1/2, and this value coincides with the phase tran- 
sition potential energy of the model for a £ [0,1]. At 
the same time, the value v = 1/2 is also the maximum 
potential energy per spin of the model, and the phase 
transition is of the "partial equivalence of ensembles" - 
type as described in |21| . For a > 1 , Af converges to the 
finite value l/(2C(a)) in the thermodynamic limit, and 
the resulting bound (|38|) is strictly positive and cannot 
give any indication of the phase transition of the model. 



VI. DISCUSSION AND CONCLUSIONS 

We have analyzed the stationary points of the poten- 
tial energy function ([5]) of the classical XY chain with 
power law pair interactions, decaying like r~ a with the 
distance r on the lattice. Computing all stationary points 
of V seems to be way out of reach, but special classes of 
stationary points can be constructed. For the class of 
"spinwave-type" stationary points where all differences 
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9i — 6{-i between neighboring angles are equal, we have 
analytically computed, in the limit of large system size 
N, the asymptotic behavior T> a of the Hessian determi- 
nant of V as a function of the potential energy per spin 
v. The computation is based on the Toeplitz property of 
the Hessian and makes use of a Szego-type theorem. The 
analytic results have been compared to numerical com- 
putations of the Hessian determinants for system sizes 
of up to N = 501, and the agreement was found to be 
excellent. 

The motivation behind these calculations is based on 
a recently discovered relation between phase transitions 
and the properties of stationary points of classical many- 
body Hamiltonian functions, as reviewed in Sec. HI Ac- 
cording to this relation, a phase transition is signaled by 
the vanishing of the (suitably scaled) Hessian determi- 
nant of V, evaluated along a suitably chosen sequence of 
stationary points of V, in the thermodynamic limit 
Moreover, the thermodynamic limit value of the (poten- 
tial) energy of such a sequence of stationary points coin- 
cides with the phase transition (potential) energy of the 
model described by V. 

For the XY chain with power law pair interactions 
with exponent ^ a ^ 1, we found that the asymp- 
totic value T> a of the Hessian determinant at the spinwave 
stationary points is zero at the potential energy per spin 
v = 1/2. In agreement with the criterion on phase transi- 
tions and stationary points, this value coincides precisely 
with the phase transition potential energy of the model. 

For a > 1, T> a is bounded away from zero, giving no in- 
dication of the phase transition occurring for 1 < a 2. 
This is of course a somewhat disappointing result, as this 
is the most interesting case: For a > 1 an exact solution 
of the thermodynamics of the XY chain is not known, 
and obtaining an exact expression for the critical en- 
ergy would have been a remarkable result. However, the 
reader should keep in mind that we have considered only 
the special class of spinwave stationary points. One good 
reason to focus on this class was the observation that, for 
the XY chain with nearest-neighbor interactions studied 
in Ref . [l7j , the spinwave stationary points were the "flat- 
test" ones (in the sense of having the smallest value of 
T> a for a given value of the potential energy per spin v), 
and therefore good candidates for V a to vanish. On the 
other hand, the presence of spinwave stationary points 
depends crucially on the boundary conditions, and their 
number grows slower than exponentially with the system 
size. Hence, in order to find an asymptotically vanishing 
Hessian determinant in the sense of Q , it appears to be 



necessary to go beyond the study of spinwave stationary 
points. 

Beyond the analysis of specific features of the XY 
chain with power law interactions, the results reported 
in the article provide a number of more general indica- 
tions that might prove useful for further applications of 
the criterion of Sec. HI First, the strategies of how to con- 
struct special classes of stationary points can be extended 
straightforwardly to other types of spin-spin-interactions 
and to higher-dimensional lattices. Second, Szego-type 
results should be applicable for the computation of the 
large- N asymptotics of the Hessian determinant also in 
other one-dimensional models. And third, the crucial 
step for successfully applying the criterion to some model 
is certainly the choice of a suitable class of stationary 
points. The present case study, although not fully con- 
clusive by itself, provides a further piece of information 
that can contribute toward an understanding of this is- 
sue. 

An interesting open question is why the spinwave sta- 
tionary points discussed in this article allow one to suc- 
cessfully detect the phase transition for exponents a be- 
tween zero and one, but fail in the case of a € (1, 2], but 
unfortunately the author can only speculate about the 
reasons. In general, one would expect that a phase tran- 
sition can be detected successfully only from stationary 
points that are somehow "relevant" for the transition, 
in the sense that they correspond to states that domi- 
nate the behavior of a phase in the thermodynamic limit. 
Physical intuition about the phases of a model should 
therefore be of help when choosing a class of stationary 
points. But for the XY model with power law interac- 
tions, no significant difference seems to distinguish the 
cases a e [0, 1] and a <E (1, 2]: In both cases, the transi- 
tion separates a ferromagnetically ordered from a para- 
magnetic phase, and it remains unclear to the author why 
spinwave stationary points capture the transition in one 
case, but not the other. 
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